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CHARACTERISTICS  OF  PLASMA-BEAM 


WALL  INTERACTIONS 


General  Remarks 

The  present  report  covers  work  performed  under  Grant  AFOSR-78-3603 
concerning  a  series  of  investigations  related  to  high-density  plasma 
beams. 

The  general  purpose  of  the  work  was  aimed  at  understanding  the 
physical  characteristics  of  the  beam  in  order  to  use  the  beam  in 
several  applications  of  great  interest  to  the  Air  Force.  Most  of  these 
required  the  formation  and  injection  of  dense,  highly  ionized  plasmas 
in  magnetic  "compression"  machines.  These  machines  then  yield 
according  to  the  extent  of  compression  used  X-rays  or  neutrons  fluxes 
or  both. 

The  production  of  hot  plasmas  of  high  density,  of  course,  has 
been  known  for  several  years  following  the  classic  work  on  the  plasma 
focus.  The  ability  to  produce  hot  plasma  beams  of  very  high  density, 
however,  is  not  readily  found,  specially  if  one  requires  the  beam  to 
be  well  collimated,  dynamically  stable  and  with  densities  better  than 
10  /cm  . 

The  work  described  here,  substantially  summarizes  the  effort  of 
several  years.  It  is  felt  that  the  state  of  the  art  in  the  design  of 
plasma  guns  of  the  Cheng  type  is  finally  understood.  We  cure  also 
convinced  that  because  of  this  understanding  one  can  construct  a 
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group  of  guns  to  be  operated  collectively  for  the  purpose  of  injecting 
plasmas  of  several  configurations  and  this  with  a  high  degree  of 
reliability. 

The  work  described  in  the  report  includes  the  contribution  of 
several  investigators  in  particular  Dr.  A.  Ferendeci  and  Mr.  R.  Mesli 
and  to  a  lesser  extent  J.  Kasky  and  R.  Webster. 

Part  of  these  investigations  have  been  reported  at  the  1979  IEEE 
Conference  on  Plasma  Science,  Montreal,  Canada  (Papers  IB  11).  A  series 
of  papers  are  presently  being  proposed  for  submission  to  various 
Scientific  Journals. 
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CHAPTER  I 


INTRODUCTION 


The  Shiva  machine  is  a  plasma  device  used  to  generate  intense 
X-rays1  by  compressing  a  plasma  shell  by  the  magnetic  self  field  of 
a  current  passing  through  the  shell. 

The  technique  originally  used  to  produce  the  plasma  shell  was 
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to  implode  a  thin  cylindrical  A£  foil  '  in  the  implosion  chamber  by 
discharging  a  high  energy  capacitor  bank.  In  this  way,  the  tempera¬ 
ture  density  and  stability  requirements  of  the  Shiva  machine  were 
fulfilled.  One  of  the  experimental  difficulties  was  that  the  implo¬ 
sion  chamber  had  to  be  replaced  after  each  shot  because  of  mechanical 
and  radiation  damage. 

To  avoid  the  difficulty  described  above,  an  alternative  techni¬ 
que  was  suggested  to  replace  the  A £  foil  by  a  cylindrical  puff  of 
4 

gas.  The  gas  was  injected  into  the  implosion  chamber  through  cir¬ 
cular  openings  in  one  of  the  electrodes.  A  cylindrical  plasma  shell 
is  then  formed  by  discharging  the  condenser  bank  across  the  chamber. 
The  concept  of  substituting  a  gas  shell  in  place  of  the  liner  foil 
offers  definitive  advantages.  In  Ref.  4,  extensive  investigations 
have  been  carried  cut  towards  the  under standing  of  the  dynamics  of 
the  puffed  gas  mode.  However,  the  experimental  results  from  this 
technique  indicated  a  non-uniform  density  distribution  of  the  gas 


shell  both  in  the  axial  and  azimuthal  directions. 


A  second  alternative  approach  to  substitute  the  liner  foil  in 
the  Shiva  implosion  chamber  by  a  plasma  shell  rather  than  an  unionized 
gas  shell  was  proposed  a  few  years  ago  by  M.  Wolfe.^  The  plasma, 
which  has  to  be  produced  initially  by  a  coaxial  gun,  is  then  injected 
into  the  implosion  chamber  through  circular  holes.  Hopefully,  a 
cylindrical  uniform  plasma  shell  will  be  formed  between  the  plates 
of  the  chamber  which  will  consequently  be  imploded  by  the  capacitor 
bank.  From  the  requirements  of  the  Shiva  machine  100  y  gm  of  plasma 
is  needed  to  be  injected  into  the  implosion  chamber.  This  is  accom¬ 
plished  through  90  holes,  each  having  a  diameter  of  1  cm  and  aligned 
around  a  circle  having  20  cm  radius.  The  corresponding  minimum  number 

densities  for  different  gases  satisfying  the  Shiva  criteria  were 
0 

found  to  be: 


Hydrogen 

17 

4.23x10 

-3 

can 

Helium 

1.06xl017 

-3 

cm 

Aluminum 

1.57xl016 

-3 

cm 

Copper 

6.66xl015 

-3 

cm 

TABLE 

1.1 

The  separation  between  the  plates  of  the  implosion  chamber  is  2  cm, 
therefore  the  uniformity  of  the  injected  plasma  in  both  the  axial 
and  azimuthal  directions  are  of  great  importance.  Also,  during  the 
injection  period  there  should  be  very  small  recombination  of  the 
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injected  plasma  such  that  the  required  mass  density  earn  be  maintained. 

In  summary,  then,  the  idea  of  the  plasma  injection  into  the 

Shiva  machine  consisted  of  producing  a  symmetric,  highly  dense  and 
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energetic  cylindrical  plasma  shell  from  a  subsidiary  gun  '  and  of 
directing  this  plasma  into  the  Shiva  implosion  chamber.  Theory  pre¬ 
dicts  that  any  initial  fluting  in  the  plasma  will  damp  out  during  the 
implosion.^ 

The  work  described  in  this  dissertation  was  for  the  purpose  of 
testing  this  concept.  A  plasma  beam  produced  by  a  small  scale  co¬ 
axial  gun,  was  to  be  injected  through  a  hole  drilled  in  the  center 
of  a  plate.  A  second  solid  plate  placed  2  cm  a  .-/ay  from  the  first  one 
simulated  the  unperforated  electrode  of  the  implosion  chamber.  In 
this  study,  special  attention  was  to  be  directed  to  the  production 
of  an  efficient  injection  scheme  and  to  the  understanding  of  the 
dynamics  of  the  plasma  flowing  through  openings  in  a  plate. 

The  incident  plasma  beam  is  produced  by  an  improved  coaxial 
g 

Cheng  gun.  Various  diagnostic  techniques  are  reported  here  later  to 
study  the  behavior  and  the  dynamics  of  the  plasma  beam  before  and 
after  passing  through  the  hole.  The  spatial  distribution  of  the  in¬ 
jected  plasma  in  both,  axial  and  azimuthal  directions  was  measured. 

The  important  parameters  that  have  been  determined  were:  the 
temperature,  density  and  directed  velocity  of  the  plasma  beam.  The 
diverging  flow  of  the  jet  as  it  comes  out  of  the  plate  has  been 
observed.  Electron  density  measurements  were  initially  made  by  means 


of  laser  interferometry.  These  measurements  were  checked  by  optical 
spectroscopy  which  also  yielded  the  temperature.  High  speed  photo¬ 
graphic  techniques  gave  velocity  measurements  and  aided  in  the  obser¬ 
vation  of  the-  plasma  behavior .  The  formation  of  the  plasma  beam  was 
monitored  by  the  time  dependent  voltage  and  current  characteristics 
of  the  discharge  mechanism. 

Concurrently,  theoretical  work  using  a  single  fluid  model  was 
undertaken  in  an  attempt  to  analyze  the  flow  pattern  of  the  jet  as 
it  escaped  into  the  implosion  chamber.  This  is  given  in  details  in 
Chapter  2. 

A  description  of  the  experimental  layout  and  the  diagnostic 
techniques  used  to  perform  the  measurements  mentioned  above  is  out¬ 
lined  in  Chapter  3. 

In  the  fourth  chapter  the  experimental  results  of  the  observa¬ 
tions  made  are  discussed. 

Finally,  in  the  fifth  chapter  the  theory  is  compared  with  the 
experimental  observations  and  suggestions  for  further  study  are  made. 


CHAPTER  II 


< 


THEORETICAL  CONSIDERATIONS 

2.1  Formulation  of  the  Problem 

The  theoretical  discussion  in  this  chapter  deals  with  the 
dynamics  of  the  plasma  injection  mechanism  into  the  Shiva  implosion 
chamber . 

To  understand  the  justification  of  the  assumptions  underlying 
the  equations  of  motion  for  the  plasma,  a  few  comments  on  the  construc¬ 
tion  of  the  implosion  chamber  are  in  place. 

The  device  consists  of  two  parallel  electrodes  (Fig.  2.1),  one 
of  which  is  perforated  to  allow  the  injection  of  the  plasma.  The 
condenser  bank  used  in  conjunction  with  the  Shiva  machine  is  allowed 
to  discharge  shortly  after  the  injection  has  taken  place.  At  that 
instant  the  plasma  is  concentrated  mainly  in  the  vicinity  of  the  entry 
holes.  The  arc  discharge  will  be  initiated  at  the  place  of  minimum 
inductance  in  the  plasma  that  is  providing  a  galvanic  conduction  path 
for  the  current. 

Therefore,  the  conditions  to  be  investigated  will  correspond  to 
the  behavior  of  the  injected  plasma  in  the  vicinity  of  the  holes. 

This  can  be  approximated  quite  reasonably  by  a  steady  flow  of  the 
plasma  as  it  passes  through  the  hole.  Indeed,  one  notices  that  the 
transit  time  of  the  plasma  through  the  thickness  of  the  plate  is  very 
short  compared  to  the  lifetime  of  the  plasma  beam  incident  on  the 


plate. 
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Plasma 

Gas  Reservior 


Second  electrode 


implosion  chamber 


perforated 

electrode 


1.  Schematic  representation  of  the  idealized 
geometry  for  the  plasma  injection. 


Because  the  arc  discharge  occurs  very  soon  after  the  incidence 
of  the  plasma  in  the  chamber  and  before  the  flow  is  appreciably  altered 
by  the  back  electrode,  we  can  neglect  the  effect  of  this  plate  on  the 
flow  pattern  of  the  jet. 

The  model  used  to  study  the  injection  mechanism  assumes  thus 
that  a  steady  state  He  plasma  beam  is  incident  on  a  plate  placed  per¬ 
pendicular  to  the  axis  of  the  beam.  The  plate  has  a  single  infinite 
slit  of  width  a  as  shown  in  Fig.  2.2.  The  slit  width  is  assumed  to 
be  smaller  than  the  diameter  of  the  beam  D.  We  assume  that  the  results 
obtained  for  a  single  slit  can  be  generalized  to  the  multiple  holes 
of  the  Shiva  machine. 

Earlier  experimental  results  on  the  plasma  beam  showed  that  the 
beam  is  radially  uniform,  very  dense  and  that  it  propagates  with 
supersonic  speeds.10'11  Observations  have  also  shown  that  the  jet 
possesses  an  irrotational  motion  and  can  be  considered  a  compressible 
fluid.  As  a  first  approximation,  a  single  fluid  model  is  chosen  for 
the  plasma  by  assuming  a  small  value  of  Debye  Length.  The  molecular 
effects  on  the  plasma,  i.e.,  recombination,  viscosity,  heat  conduction 
and  heat  transfer  are  neglected.  Formation  of  shocks  are  excluded  in 
this  model.  Magnetic  forces  are  excluded  by  assuming  that  there  are 
no  trapped  fields  in  the  plasma  and  that  no  currents  flow  between  the 
gun  and  the  perforated  plate. 

The  validity  of  these  assumptions  have  been  verified  by  the 
experimental  results  obtained  through  the  measurements. 


A  two  dimensional  model  is  considered  as  shown  in  Fig.  2.2.  The 


x-axis  is  taken  to  be  along  the  axis  of  the  flow  and  the  y-axis  per¬ 
pendicular  to  it.  The  width  of  the  slit  in  the  z  direction  is  assumed 
infinite. 

Most  of  the  work  found  in  the  literature,  on  the  theory  of  super 
sonic  gas  flow  from  an  orifice  has  dealt  with  neutral  compres'  ble 
gases  flowing  through  holes  and  nozzles.  In  most  of  the  cases  dis¬ 
cussed,  the  back  pressure  is  comparable  to  that  of  the  incident  pres- 
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sure  of  the  gas.  '  In  the  present  work,  the  ratio  of  the  exit 
to  the  back  pressure  of  the  outer  medium  is  very  high  since 
the  jet  flows  into  a  high  vacuum.  The  plasma  flow  is  still  expected 
to  form  a  diverging  jet  structure. 

2.2  Governing  Equations 

Using  the  above  assumptions,  the  governing  equations  for  the 

escaping  jet  are  deduced  from  the  general  magnetohydrodynamic  equa- 

. .  .  14,15 

tions  and  are 


p(V*V)V  +  7p  = 

0 

(2.1) 

V* (pV)  = 

0 

(2.2) 

=  y  E 
dp  p 

c2 

(2.3) 

Vxv 

0 

(2.4) 

ds  =  C  ^  - 

c 

dp 

(2.5) 

V  p 

p 

P 

where  p  =  mass  density  =  nm 

•> 

V  =  fluid  velocity  vector  =  average  flow  velocity 


p 


kinetic  pressure  =  n  k  T 


Y  =  ratio  of  specific  heats  =  C  /C 

P  V 

1/2 

C  =  speed  of  sound  =  (y  kT/m) 
s  =  entropy 


The  five  equations  (2.1)  to  (2.5)  express  respectively  conservation  of 
momentum,  conservation  of  mass,  the  equation  of  state,  the  condition 
for  irrotational  flow  and  the  equation  of  state  for  an  ideal  gas  flow. 

Combining  equations  (2.1),  (2.3)  and  (2.5)  into  equation  (2.2), 
the  governing  equations  of  a  supersonic  steady  continuous  homoenergetic 
and  irrotational  flow  of  an  ideal  gas  with  constant  specific  heats 
reduce  to1"* 

2 

-y  1  -+  V 

v-v  =  ~  V-V  (  -  )  (2.6) 

C 

Vxv  =  0  (2.4) 

For  a  cartesian  system,  equations  (2.4)  and  (2.6)  become 


2  2 
u  .  3u  uv  3u  uv  5u  v  3v 

c2)  3x  "  c2  3y  ~  c2  3x  +  U  “  Q2]  3y 


0 


(2.7) 


3v  3u 

3x  3y 


(2.8) 


where  u  and  v  denote  the  components  of  the  flow  velocity  V  in  the  x 
and  y  directions,  respectively. 


2. 3  Method  of  Characteristics 

Expressions  (2.7)  and  (2.8)  form  a  set  of  non-linear  partial 
differential  equations  of  the  hyperbolic  type.  The  system  of  equations 


/ 
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is  solved  by  the  well  known  numerical  method  of  characteristics.  Ref. 
15,16,17,18,19  and  20  deal  extensively  with  the  method  of  character¬ 
istics. 

The  solutions  of  equations  (2.7)  and  (2.8)  represent  the  slope 
of  the  physical  characteristics  (Mach  waves)  at  each  point  of  the 
flow  and  are  given  in  Ref.  15  as: 


Making  use  of  the  relationship  between  the  magnitude  of  the  velocity 
V  and  its  components  u  and  v  and  the  flow  angle  8  (see  Fig.  2.3)  as 
well  as  the  relationship  between  the  Mach  angle  a  and  the  Mach 
number  M  (see  Fig.  2.4),  an  alternate  form  of  expressions  (2.9)  and 
(2.10)  can  be  written  as 


C 


+  ' 


*2.  m 

dx 


tan  (9  +  a) 


C  : 


§2.  = 
dx 


(2.11) 


tan  (8  -  a) 


(2.12) 


TABLE 


/ 
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Relationship  between  the  velocity  com¬ 
ponents  u  and  i.  and  the  velocity  magnitude  V'  and  the 
flow  angle  8  (Ref.  17) 


Relationship  between  the  Mach  angle  > 
and  the  Mach  number  M  (Ref.  17) 


Figure  2.3  Figure  2.4 


y 

A 
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Here  the  flow  angle  6  is  defined  as 


9  =  tan  (— ) 

v 


(2.13) 


and  the  Mach  single  a  as 


-1  C  -1  1 

“  =  Sin  -  -  Sin  - 


(2.14) 


where  the  Mach  number  M  is  defined  in  the  usual  manner  as 


M 


V 

C 


(2.15) 


The  characteristics  C+  and  C  as  the  Mach  lines  at  a  point  in  a  flow 
field  are  shown  in  Fig.  2.5. 

The  equations  relating  the  changes  in  direction  of  flow  measured 
by  0  and  the  resulting  change  in  the  magnitude  of  the  gas  velocity  V 
across  Mach  waves  provide  a  workable  expression  for  the  numerical 
method  of  characteristics.  These  can  be  written  as*^'^  : 


d9  =  + 


J7~ 

V  M  -  1 


M(1  +  M2) 


dJM 


(2.16) 


The  positive  and  negative  forms  of  equation  (2.16)  refer  to  the 
crossing  of  a  Mach  wave  whose  slopes  are  (dy/dx)  =  tan  (9  -  a)  and 
(dy/dx)  =  tan  (0  +  a)  ,  respectively.  The  familiar  form  of  equation 
(2.16)  is 

de  = 


+  dv 


(2.17) 


For  a  perfect  gas,  the  integrated  form  of  (2.17)  is  written  as: 


21 

v (M)  is  called  the  Prandtl  Meyer  angle  and  is  defined  as  the 
angle  through  which  the  stream  flow  turns  in  expanding  from  a  Mach 
number  of  1  to  a  supersonic  Mach  number  M.  Equation  2.17  is  the  main 
expression  for  the  step  by  step  numerical  integration  procedure  used 
in  the  method  of  characteristics. 

A  complete  set  of  tables,  equations,  and  figures  to  be  used  in 
the  analysis  of  Helium  gas  flowing  at  supersonic  and  hypersonic 
speeds  are  given  in  Ref.  21. 


2.4  Application  for  a  Flow  with  M  =  2 

The  flow  in  the  vicinity  of  the  opening  is  obtained  as  stated 
in  section  2.1,  by  neglecting  the  presence  of  the  second  plate.  From 
the  experimental  results,  the  incident  plasma  beam,  before  it  reaches 
the  perforated  plate  has  a  Mach  number  approximately  2.  Therefore, 
the  method  of  characteristics  have  been  applied  for  the  case  of  an 
initial  Mach  number  M  =2.  From  symmetry  considerations,  calcula¬ 


tions  have  been  carried  out  only  for  the  upper  half  of  the  flow  field. 


The  corresponding  computations  are  given  in 

Table  2.1.  The  notation  04  used  is  defined^-5  as  the  flow  angle 
corresponding  to  M  =  1.  Since  for  M  =  1  the  value  of  v  becomes  zero. 


the  integration  of  equation  (2.15)  gives  the  following  expressions: 
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6  -  =  +  v (M) 

where  the  superscripts  (-)  and  (+)  refer  to  waves  whose  Mach  waves 
have  slopes  dy/dy  =  tan(0  +  a),  respectively. 

The  underlined  numbers  refer  to  the  quantities  known  at  the 
beginning  of  the  computation  of  the  other  quantities  in  each  row. 
Throughout  all  the  computations,  Helium  is  assumed  to  be  a  perfect 
plasma  gas. 

The  Mach  wave  pattern  constructed  are  illustrated  in  Fig.  2.6. 

Figure  2.7  shows  the  plot  of  the  density  and  temperature  variations 

of  the  escaping  jet  in  the  x-direction  for  two  vertical  distances 

(y  =  0  cm  and  0.25  cm)  away  freon  the  axis  of  symmetry.  In  Fig.  2.8, 

the  variation  of  the  same  parameters  in  the  vertical  direction  are 

plotted  for  the  planes  x  =  1  cm  and  x  =  2  cm  away  from  the  orifice. 

*  * 

In  Figures  2.7  and  2.8,  the  dimensionless  ratios  P/Pt  and  T/T^  have 
been  utilized  rather  than  the  absolute  density  p  and  temperature  T 
in  order  to  make  the  results  more  general.  In  this  way,  the  results 
are  independent  of  the  initial  beam  density  and  temperature. 

The  theoretical  model  predicts  that  at  the  vicinity  of  the 
hole,  the  injected  plasma  flow  maintains  the  same  density  and  tempera¬ 
ture  of  the  incident  plasma.  It  is  found  that  the  initial  state 
extends  within  a  triangular  region  formed  by  the  initial  Mach  lines 
originated  from  the  two  corners  of  the  opening.  As  the  injected  flow 


The  subscript  t  refers  to  the  total  conditions  that  would  exist  if 
the  gas  were  brought  to  rest  isentropically^ . 


for 
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expands  further  down  the  orifice,  the  density  and  temperature  of  the 
flow  decrease  significantly,  but  it  is  found  that  within  2  cm  away 
from  the  hole,  the  total  mass  density  of  gas  required  for  the  Shiva 
experiment  is  fulfilled. 

2. 5  General  Behavior  with  Mach  Number  M 

Due  to  the  value  of  the  initial  Mach  number  M,  the  Mach  wave 
pattern  for  the  jet  has  to  be  modified.  A  plot  predicting  the 
dependence  of  the  initial  Mach  wave  angle  6^  and  the  escaping  flow 
angle  02  on  the  Mach  number  M  is  represented  in  Figure  2.9.  The 
change  in  the  flow  angles  0^  characterizing  the  same  parameters  of 
the  initial  state  of  the  plasma  beam  affects  the  spatial  extent 
of  the  triangular  region  as  '•.hown  in  Fig.  2.10.  It  is  concluded 
that  high  Mach  numbers  will  produce  smaller  initial  angles  extending 
the  initial  plasma  state  further  down  the  hole.  Also,  the  change  in 
the  escaping  flow  angle  8 2  as  the  initial  Mach  number  M  is  varied 
causes  a  change  in  the  extent  X2  of  the  plasma  mixing  region  when 
more  than  one  hole  is  considered  in  the  plate.  Fig.  2.11  shows  a 
plot  of  the  dependence  of  the  extent  X2  on  the  initial  Mach  number  M. 
One  notices  that  for  low  Mach  numbers,  due  to  the  presence  of  multiple 
holes,  the  plasma  mixing  takes  place  at  the  vicinity  of  the  holes 
while,  at  higher  Mach  numbers,  no  mixing  region  from  the  neighboring 
holes  is  expected. 

Optimizations  have  to  be  made  in  order  to  get  both  large  extent 


of  plasma  initial  state  and  plasma  mixing  region  in  the  Shiva 


MACH  NUMBER  M 


Figure  2.10.  Continuation  of  the  plasma  initial  state  on 


the  Mach  Number  M. 
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implosion  chamber. 

A  more  appropriate  solution  to  the  proposed  injection  problem  is 
to  use  a  circular  hole  in  the  plate  instead  of  one  dimensional  slit 
and  use  cylindrical  coordinates  for  the  solution  of  equations  (2.4) 
and  (2.6) ,  but  the  two  characteristics  on  the  hodograph  plane  and  the 
physical  plane  are  not  decoupled,  as  for  the  rectangular  case15  and 
thus  requires  computer  solutions  to  perform  the  many  iterations 
necessary  to  solve  the  problem.  A  computer  program  is  not  presently 
available  to  carry  out  these  calculations  and  should  be  the  subject 
of  further  study. 


) 
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CHAPTER  III 

EXPERIMENTAL  ARRANGEMENT  AND  DIAGNOSTIC  TECHNIQUES 

The  experiment  and  the  diagnostic  techniques  performed  were 
chosen  to  provide  a  means  of  testing  the  feasibility  of  the  Shiva 
plasma  injection  concept. 

The  plasma  beam  generated  as  part  of  this  work  was  also  inves¬ 
tigated  and  characterized  to  determine  its  appropriateness  to  testing 
the  Shiva  plasma  injection  concept. 

3 . 1  Experimental  Arrangement 

3.1.1  General  description 

The  beam  is  produced  with  an  improved  version  of  Cheng's  plasma 
9 

gun,  a  coaxial  gun  of  conical  shape  with  0.7  cm  inner  electrode 
radius  and  ^  4  cm  mean  outer  electrode  radius.  The  gun  was  contained 
in  a  Pyrex  glass  tube  where  most  of  the  experiments  were  carried  out. 
The  whole  system  was  evacuated  to  a  pressure  of  ^  10  6  Torr. 

The  rest  of  the  experimental  apparatus  was  comprised  of  the  main 
energy  storage  capacitor  bank,  an  electromagnetic  valve  for  gas 
puffing  with  related  energy  storate  capacitor,  the  charging  high  vol¬ 
tage  system,  the  electronic  triggering  units  and  various  diagnostic 
equipment.  An  illustration  of  the  experimental  arrangement  is  shown 
in  Figure  3.1.  The  schematics  of  the  energy  storage  bank  and  control 
system  for  the  plasma  gun  as  well  as  the  schematic  of  the  trigger 
unit  can  be  found  in  Ref.  10. 
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3.1.2  Modifications 


The  gun  is  operated  in  the  puff  mode.  In  order  to  achieve  the 
minimum  densities  for  the  Shiva  experiment,  estimated  to  be  'v  1017  cm" 
for  the  He  gas  (Table  1.1),  the  original  electromagnetic  puff  valve10,1'* 
had  to  be  redesigned.  The  new  version  allowed  more  gas  to  be  puffed 
into  the  gun. 

While  the  control  system  remained  the  same  as  the  one  used  by 
Lee10,  the  energy  storage  system  was  revised  for  a  better  match  of  the 
capacitors  to  the  gun  thus  assuring  more  efficient  energy  transfer. 

In  particular,  the  gun  switch  has  been  removed  altogether  to  reduce 
the  system  inductance. 

A  conducting  back  plate  has  been  added  to  the  original  version 
22  23 

of  Cheng's  gun.  '  This  configuration  causes  faster  breakdown  of 
the  gun  and  allows  lower  voltage  operation.  It  also  improved  the 
reproducibility  of  the  beam  formation  considerably.  A  schematic  of 
the  coaxial  gun  presently  used  is  given  in  Figure  3.2. 

In  addition,  many  improvements  have  been  made  on  electrical 
protection  of  the  gun  and  the  valve.  Particular  attention  was  given 
to  grounding  and  shielding  problems.  The  negative  polarity  of  the 
inner  electrode  has  been  determined  experimentally  to  be  fundamental 
to  the  beam  formation. 

3.1.3  Mode  of  operation 

The  procedure  of  firing  the  gun  is  as  follows:  a  15yF  energy 
storage  capacitor  is  used  to  control  the  fast  electromagnetic  valve. 


Figure  3.2.  Schematic  of  the  Coaxial  Plasma  Accelerator. 
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The  back  pressure  of  the  valve  and  the  on-time  of  the  electromagnetic 
valve  determine  the  amount  of  He  gas  puffed  into  the  barrel  of  the 
gun.  A  60yF  low  inductance  capacitor  bank  charged  to  the  proper  vol¬ 
tage  was  permanently  connected  across  the  two  electrodes  of  the  gun. 

An  initial  manually  triggered  pulse  starts  the  valve  capacitor 
to  discharge  causing  a  momentary  opening  of  the  valve  allowing  He  gas 
to  enter  the  barrel  of  the  gun.  The  breakdown  of  the  gun  occurred 
when  the  Paschen  conditions  are  reached.  The  time  elapsed  between  the 
firing  of  the  puff  valve  and  the  triggering  of  the  gun  is  found  to 
depend  on  the  capacitor  voltages  and  the  background  pressure  of  the 
valve . 

The  mechanism  of  the  beam  formation  has  been  explained  as  a 

J*B  force  acting  upon  the  electrons,  while  the  ions  experienced  the 

22 

Coulomb  force  acceleration  due  to  the  ejected  electrons.  The  cur¬ 
rent  sheet  is  therefore  accelerated  down  to  the  muzzle  of  the  gun  and 
then  collapses  in  front  of  the  inner  electrode  resulting  in  the 
ejection  of  a  high  density  high  velocity  plasma  beam. 

3.1.4  Simulation  of  the  injection  experiment 

Initially,  an  interaction  chamber  made  of  brass  material  was 
inserted  in  the  experimental  set-up  in  order  to  study  the  beam-plate 
interaction.  The  presence  of  the  interaction  chamber,  its  shape  and 
grounding  have  considerably  changed  the  properties  of  the  beam,  pro¬ 
ducing  appreciable  return  currents  and  modifying  the  beam  charac¬ 


teristics. 


In  the  next  stage,  the  interaction  chamber  was  removed  and  re¬ 


placed  with  Pyrex  glass  tubing. 

For  the  interaction  plates,  initially  A1  or  Cu  plates  were  used. 
It  was  found  out  that  the  energetic  plasma  beam  ablates  a  considerable 
amount  of  material  from  those  electrodes.  The  light  generated  from 
these  ablatants  and  the  line  impurities  made  the  study  of  the  beam 
hole  interaction  almost  impossible.  Therefore,  these  plates  were  re¬ 
placed  by  two  tungsten  plates  separated  by  0.5  cm*  simulating  a  one 
dimensional  slit  representing  the  actual  assumed  theoretical  model. 

3.2  Diagnostic  Techniques 

The  diagnostic  procedures  were  directed  towards  the  analysis  of 
the  beam  features  before  and  after  the  opening. 

The  diagnostic  techniques  performed  to  infer  the  different 
plasma  parameters  were  mainly  optical  interferometry  optical  spectro¬ 
scopy,  high  speed  visual  photography  and  monitoring  of  the  electrical 
parameters  of  the  gun  system.  These  will  be  described  separately 
below. 

3.2.1  Optical  interferometry 

Within  the  range  of  the  electron  densities  produced  by  the  gun 
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an  Ashby-Jephcott  type  laser  interferometer  was  found  to  be  suitable 
for  density  measurements.  The  visible  radiation  of  a  He-Ne  laser  at 

O 

6328A  was  used  in  the  interferometer. 


3.2.2  High  speed  photography 


High  speed  visual  photography  measurements  were  made  using  an 
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STL  image  converter  camera  both  in  the  framing  and  the  streak  mode. 

The  camera  is  triggered  either  optically  or  electrically. 

3.2.3  Optical  spectroscopy 

The  spectroscopic  techniques  used  were  concerned  with  the  esti- 

O  O 

mation  of  the  Hei(5876A)  and  He^ (4686A)  line  intensities  as  well  as 
the  continuum  intensity. 

These  measurements  were  performed  by  a  McPherson  (Model  2165) 

spectrometer.  The  integrated  line  spectra  is  recorded  on  a  polaroid 

film.  Slits  to  monitor  the  time  history  of  the  spectra  corresponding 

to  He^,  He^  and  continuum  were  machined  on  phosphor  bronze  metal 

strip  and  placed  at  the  plane  of  the  exit  slit.  Three  photomultipliers 

were  placed  behind  each  slit  and  care  is  taken  to  match  the  outputs 

of  the  photomultiplier  tubes  to  the  characteristic  impedances 

of  the  coaxial  lines.  The  signals  were  monitored  by  a  Tektronix  551 

dual  beam  oscilloscope.  A  reference  He  gas  discharge  source  is  used 

for  locating  the  various  lines.  The  apparatus  was  calibrated  by 
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means  of  a  carbon  arc  standard 

Since  the  time  delay  between  the  initial  triggering  of  the  valve 
and  the  breakdown  of  the  gun  was  not  reproducible,  other  means  of 
initiating  the  diagnostic  equipment  such  as  the  optical  fiber  or  the 
gun  current,  were  used.  The  time  t  =  0  was  usually  obtained  by  getting 
an  optical  trigger  from  the  initial  light  output  of  the  gun,  while 
the  initial  rise  in  gun  current  was  used  to  trigger  the  oscilloscope 
for  some  of  the  measurements. 
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The  terminal  currents  and  voltages  of  the  gun  were  monitored  by 
a  Rogowski  coil  and  capacitive  voltage  divider,  respectively. 


L 
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CHAPTER  IV 

EXPERIMENTAL  RESULTS 

4.1  Plasma  Parameters 

Experimental  measurements  were  performed  to  provide  a  check  on 
the  theoretical  calculations  of  Chapter  2. 

One  of  the  important  parameters  to  be  experimentally  determined 
was  the  Mach  number  M.  As  pointed  out  in  Chapter  2,  M  is  defined  as 


M  -  V  _  V 

M  -  c  - 


^(■ykT/m) 


(4.1) 


In  the  plasma,  the  flow  velocity  V  is  defined  as: 


V  = 


p .  V.  +  p  V 

ii  e  e 


P.  +  P 
l  e 


(4.2) 


where  p  and  v  denote  the  density  and  velocity  respectively,  and  the 

subscripts  i  and  e  refer  to  ions  and  electrons.  Since  the  electron 

mass  is  much  smaller  than  the  ion  mass,  p  <<  p .  ,  then  the  beam 

e  i 

velocity  can  be  taken  equal  to  the  ion  velocity 

V  %  v.  (4.3) 

l 

Thus  the  momentum  of  the  beam  is  carried  by  the  ions.  Using  a  simi¬ 
lar  argument,  the  temperature  T  used  in  expression  (4.1)  is  described 
as  follows.  In  a  fully  ionized  ideal  gas,  the  pressure  P  of  the  gas 
is  the  sum  of  the  partial  pressures  P^^  due  to  the  ions  and  of  Pg 
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contributed  by  the  electrons. 

P  =  P.  +  P  =  nk(T.  +  T  )  (4.4) 

i  e  i  e 
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Previous  measurements  on  the  ion  temperature  have  shown  that 

the  ions  in  the  beam  are  cold,  i.e.,  T,  <<  T  ,  thus,  the  temperature 

l  e 

T  in  the  speed  of  sound  (eq.  4.1)  is  taken  to  be  equal  to  the  electron 
temperature. 

To  evaluate  the  initial  Mach  number  of  the  beam,  the  parameters 
therefore  required  to  be  determined  a  priori,  are  the  ion  velocity  and 
the  electron  temperature.  The  electron  density  is  also  of  primary 
interest  as  well,  although  it  is  not  needed  for  flow  analysis,  but 
it  is  necessary  to  verify  the  various  assumptions  made  in  the  theore¬ 
tical  consideration  as  well  as  in  spectroscopic  analysis. 

All  these  parameters  are  also  required  after  the  passage  of  the 
plasma  through  the  central  hole  in  the  plate.  The  extent  of  the 
divergence  of  the  plasma  jet  as  it  leaves  the  plate  is  also  needed. 

4.2  Beam  Structure  Measurements 

Initial  investigations  of  the  beam  parameters  are  focused  on 
the  evaluation  of  the  beam  parameters  without  the  perforated  plate. 

4.2.1  High  speed  photography 

A  typical  set  of  image  converter  camera  framing  pictures  of 
the  beam  taken  for  different  time  intervals  at  50  nS  exposure  time 
are  shown  in  Figure  4.1.  These  pictures  reveal  that  the  beam  lasts 
for  5  'v  6  ps  and  stays  uniformly  collimated  over  'v  50  cm.  The 
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diameter  of  the  beam  is  found  to  be  2  ^  3  cm. 

The  beam  velocities  are  determined  from  the  streak  pictures  using 
an  STL  camera.  A  slit  1  mm  wide  was  placed  along  the  axial  length  of 
the  pyrex  pipe  with  strips  equally  spaced  every  5  cm  for  spatial  iden¬ 
tification.  A  total  time  exposure  of  10  ps  was  used  for  the  streak 
pictures.  A  typical  set  of  streak  pictures  taken  for  different  gun 
voltages  at  a  given  back  pressure  (2  psi)  is  shown  in  Figure  4.2. 

The  valve  voltage  was  kept  constant  at  8  kv.  It  was  found  that  the 
initial  beam  front  propagates  with  a  velocity  of  10  'v  20  cm/ps  de¬ 
pending  on  the  gun  voltage  and  the  back  pressure  of  the  valve  and 
travels  v  50  cm  down  the  Pyrex  tube  verifying  the  framing  picture 
observations.  The  dependence  of  the  initial  velocity  of  the  beam  on 
the  gun  voltage  and  the  valve  back  pressure  is  plotted  in  Figure 

4.3.  It  was  found  that  at  later  stages  of  the  beam,  the  beam  velocity 
decreases.  The  velocity-time  history  of  the  beam  for  different  gun 
voltages  and  different  valve  back  pressures  is  illustrated  in  Fig. 

4.4.  Figure  4.5  represents  the  velocity-time  history  of  the  beam  for 
different  gun  voltages  but  for  a  given  valve  back  pressure  (v  1  psi) . 

In  figures  4.4  and  4.5,  the  velocities  are  normalized  to  the  initial 
velocity  V  .  The  time  dependence  of  the  beam  indicates  that  the 
velocity  decreases  constantly  towards  the  trailing  edge.  The  streak 
picture  also  verify  that  the  lifetime  of  the  beam  is  v  5  ps  consistent 
with  the  framing  pictures. 

At  low  gun  voltages  the  beam  parameters  were  not  easily  repro¬ 


ducible  indicating  an  erratic  form  for  the  beam. 
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Figure  4.5.  Velocity  time  variation  of  the  beam 


given  back  pressure  1  psi. 
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For  later  analysis  on  the  understanding  the  dynamics  of  the  beam 
formation  terminal  current  and  voltage  characteristics  of  the  gun  were 
also  recorded.  Typical  current  and  voltage  waveform  for  the  gun  are 
shown  in  Figure  4.6. 

4.2.2  Spectroscopic  techniques 
a)  General  Remarks 

The  literature  deals  extensively  with  both  the  theory  and 
experiments,  on  the  radiated  intensities  of  the  He^  He^  lines  as 
well  as  the  continuum  intensity.  Theoretical  curves  relating  various 
ratios  of  these  intensities  as  a  function  of  temperature  are  given 
in  Ref.  26. 

The  spectroscopic  techniques  used  to  determine  the  plasma  para¬ 
meters  are  based  upon  the  assumption  of  local  thermodynamic  equili- 
26  27  28  29 

brium  (LTE)  '  '  '  for  the  plasma.  For  this  assumption  to  be 

valid,  the  required  electron  density  should  be  higher  than  the  mini- 

26 

mum  electron  density  Ne  given  by  Griem 
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To  locate  a  region  for  continuum  intensity  measurements  free 
of  impurities,  integrated  line  spectra  of  the  beam  radiation  is  re¬ 
corded  through  a  McPherson  spectrometer  using  a  Polaroid  film.  The 

O  O 

integrated  line  spectra  of  H6JIS876A)  and  He^I(4686A)  were  also 
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4.6  Typical  current  and  voltage  waveforms  for  the  gun. 
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recorded  on  the  same  film.  The  radiation  spectrum  of  the  beam  taken 
at  10  cm  from  the  muzzle  of  the  gun  is  shown  in  Figure  4.7b.  The 

O 

pictures  revealed  that  no  Hej  (5876A)  line  was  recorded  even  when  25 

shots  were  superimposed  on  the  same  film.  This  suggested  that  the 

plasma  beam  had  a  higher  electron  temperature  than  previously  re- 
23 

ported  .  For  the  continuum  a  region  free  of  impurity  lines  was 

© 

selected,  centered  at  5559A, 

b)  Temperature  Measurement 

Typical  photomultiplier  signals  giving  the  intensities  of 

O 

the  continuum  radiation  and  He.^ (4686A)  line  are  shown  in  Figure  4.8a. 
The  He^  line  intensity  was  found  to  be  of  the  same  magnitude  as  the 
corresponding  continuum  intensity  at  the  He^  location  which  was 
determined  from  a  calibrated  mask  using  a  standard  He  gas  discharge 
tube.  The  slits  in  the  mask  were  made  at  the  He^.,  He^  and  continuum 

o  o  o 

with  openings  of  14. 1A,  5.6A,  and  12. 7A,  respectively. 

The  ratio  of  the  total  intensity  of  the  He^  line  to  the  inten- 

O  O 

sity  of  the  underlying  continuum  (100A  wide)  centered  at  4686A  has 

been  employed  to  evaluate  the  beam  temperature.  From  the  curve  taken 

from  Ref.  26  and  shown  in  Figure  4.9  gives  Tg  (42  +  5  eV) .  An 

interpolating  formula  has  been  used  to  calculate  the  continuum  centered 

°  26 

at  line  from  the  5559A  continuum. 

The  time  variation  of  the  beam  electron  temperature  for  gun 
voltages  of  9  kV  and  10  kV  is  shown  in  Figure  4.10.  Higher  electron 
temperatures  have  been  recorded  for  lower  voltages.  It  was  also 
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Figure  4.8  Typical  photomultiplier  signals. 


Ratio  of  the  total  tin*  and  continuum  intaniitiai  (in  100- A  bonds  centered  at  lb*  lines) 
at  a  function  of  temperature  for  Ionlied-helium  lines.  (Kef  .  26) 
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observed  that  the  maximum  temperature  is  maintained  for  about  a  few 
microseconds. 

c)  Density  measurement 

The  electron  density  was  obtained  from  the  absolute  con¬ 
tinuum  intensity  measurements.3*3'31  A  carbon  arc  standard  was  used 
for  calibrating  the  photomultipliers. 

The  working  expression  relating  the  voltage  ratio  between  the 
pressure  of  the  carbon  arc  and  the  continuum  intensity  is  given  by 
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1.07  x  10 
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cm 
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(4.5) 


where  V  and  V  are  respectively  the  photomultipliers  output  voltages 

c  s 

o 

for  the  continuum  intensity  (5559A)  and  carbon  arc  standard. 

The  electron  density  has  been  found  to  increase  linearly 

17  -3 

with  the  gun  voltage  and  reaching  (2+0.25)xl0  cm  for  a  gun  vol¬ 
tage  of  llkV  as  shown  in  Figure  4.11.  The  density  time-history  of 
the  beam  for  different  gun  voltages  is  also  plotted  in  Figure  4.12. 

The  value  of  the  electron  density  thus  measured,  fulfills  the  require¬ 
ments  of  both  the  Shiva  experiment  (Chapter  1)  and  the  LTE  assump¬ 
tions. 


4.2.3  Optical  interferometry 

The  laser  interferometry  technique  is  initially  used  in  the 
presence  of  the  interaction  chamber  to  determine  the  electron  density 
of  the  beam.  The  relationship  between  the  number  of  fringes  N  de- 

O 

tected  and  the  electron  density  Ne  of  the  plasma  beam  for  6328A  laser 


of  the  maximum  beam  electron  density  on  the 
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radiation  is  given  by 

—  18 

N  =  5.6  x  10  N  (4.6) 

e 

Using  the  Ashby-Jephcott  interferometer,  the  following  interference 

fringes  are  obtained  as  shown  in  Figure  4.13.  From  these  interference 
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fringes,  an  average  electron  density  of  Ng  ^  4x10  cm  is  obtained. 
This  value  is  within  the  range  of  the  electron  density  values  obtained 
using  the  optical  spectroscopic  techniques. 

4. 3  Beam-Hole  Interaction 

An  Al  or  a  Cu  plate,  6.25  mm  thick,  having  a  central  circular 
hole  one  cm  in  diameter  was  placed  10  cm  in  front  of  the  gun. 

The  usual  observation  of  the  appearance  of  the  plates,  Figures 
4.14  a,b  indicated  that  considerable  electrode  material  was  ablated 
from  the  hole  region.  The  light  produced  by  the  Al  and  Cu  metal 

O  O 

vapor  had  lines  very  close  to  the  6328A  laser  radiation,  (Al  ;6336Aj 

o  o  o 

Cu:6326A;  Cu^:6325A;  Cu^ ^ : 6318A) ,  thus  interfering  with  and  making 

it  impossible  to  use  the  laser  interferometry  technique  with  the 
existing  equipment.  It  was  also  found  that  a  metal  impurity  line 
appeared  at  the  continuum  slit  making  it  impractical  to  use  the 
existing  continuum  location  for  spectroscopic  diagnostics. 

Observations  showed  that  the  ablation  from  the  interaction  plate 
was  of  the  same  order  of  magnitude  whether  the  plate  is  placed  inside 
the  interaction  chamber  or  in  a  Pyrex  tube  alone.  Therefore,  due  to 
the  possibility  of  current  loops,  the  interaction  chamber  was  removed 


1 

1 

A 


|«g  1  usec/cm 


Interferometer  output  for  the  plasma  beam 
a)  with  the  plasma 


b)  without  the  plasma. 


53 


and  only  a  Pyrex  pipe  and  an  interaction  plate  were  used. 

Later,  similar  experiments  were  repeated  for  which  the  hole  was 

replaced  by  a  5mm  slit  made  from  two  tungsten  plates  2mm  thick  as 
illustrated  in  Figure  4.14c. 

As  expected  the  presence  of  an  obstacle  in  front  of  the  beam 
reduced  the  beam  velocity.  The  velocity  of  the  escaping  jet  was  found 
to  be  approximately  5-10  cm/ys.  This  was  obtained  from  streak  pictures 
of  the  beam  after  it  emerged  from  the  slit.  Typical  streak  pictures 
shown  in  Figure  4.15  indicated  that  the  jet  traveled  a  length  of  10cm 
and  lasted  for  about  5  ys  verifying  the  steady  state  assumption 
stated  earlier.  Figure  4.16  gives  the  velocity  time-history  of  the 
beam  beyond  the  hole. 

An  attempt  was  made  to  determine  the  incident  velocity  at  the 
slit  but  due  to  the  strong  ablation  of  the  material,  and  accompanying 
intense  light,  the  STL  camera  was  always  overexposed  making  measure¬ 
ments  very  difficult. 

For  the  Cu  plate,  the  He^  line  was  observed  in  the  time  inte¬ 
grated  spectra  for  the  beam  before  the  hole  as  shown  in 

O 

Figure  4.7c.  It  was  found  also  that  a  Cu^  line  (5555A)  fell  within 
the  slit  opening  of  the  continuum,  making  temperature  and  density 
measurements  from  the  line-continuum  ratio  unreliable.  Instead, 

0  O 

the  intensity  ratio  of  the  line  He^  (4686A)  and  HeI(5876A)  was  used 
for  the  temperature  measurements.  The  dependence  of  this  ratio  on 
the  temperature  is  plotted  in  Figure  4.17.  Typical  photomultiplier 


Figure  4.15.  Typical  streak  pictures  of  the  beam 
after  the  hole. 

Total  streak  time  =  10  usee 

a)  3  psec  delay  b)  4  psec  delay 

V„  =  9  kV  V  =  8  kV  P  =  2  psi 
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Figure  4.17. 
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signals  for  the  Hej  and  He^  lines  before  and  after  the  hole  are  given 
in  Fig.  4.8b,c.  Calculations  indicated  that  the  electron  temperature 
of  the  beam  was  lowered  to  4+^  0.4  eV  by  the  presence  of  the  plate. 
Figure  4.18  gives  a  plot  of  the  electron  temperature  of  the  plasma 
jet  at  different  distances  away  from  the  orifice  and  at  different 
times,  indicating  the  consistency  of  the  measurements. 

Although  the  incident  velocity  was  not  evaluated,  it  is  reason¬ 
able  to  expect  that  since  the  temperature  remained  constant  during 
the  beam-hole  interaction,  so  would  the  velocity  of  the  beam.  The 
computation  of  the  initial  Mach  number  M  thus  became  straightforward 
and  gave  M  =  4  +_  0.4. 

The  framing  pictures  (Figure  4.19)  showing  the  beam-hole  inter¬ 
action  revealed  that  the  plasma  stagnates  as  it  impinges  upon  the 
plate  and  escapes  with  a  definite  expanding  angle  rather  than  filling 
in  the  volume  behind  the  plate.  The  stagnation  of  the  plasma  pro¬ 
vides  further  verification  of  the  fluid  approximation,  while  the 
diverging  structure  of  the  jet  verifies  the  compressibility  assump¬ 
tion  of  the  fluid.  In  fact,  the  framing  picture  (Figure  4.20)  of  the 
issuing  jet  taken  during  the  earliest  development  of  the  escape  gives 
the  escaping  flow  angle  0^  34°. 

The  framing  pictures  cited  above  were  taken  with  the  tungsten 


slit. 

For  M  =  (4  0.4),  the  theoretical  curve  of  Figure  2.8  gives 

12  <  ex  <  15°  and  35°  £  02  <  43°. 


the  plate 


Figure  4.20  Framing  picture  of  the  escaping  jet  taken 
at  2  psec.  (50  nsec  exposure) .  The  plate 
is  to  the  left. 
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CHAPTER  V 


CONCLUSION  AND  SUGGESTIONS  FOR  FURTHER  INVESTIGATIONS 

5. 1  Summary  of  Observations 

The  salient  features  of  the  beam-plasma  interaction  study  can 
be  summarized  as  follows. 

1)  The  assumption  of  a  uniformly  collimated  beam  incident 
on  the  hole  is  corraborated  by  the  framing  pictures  (Figure  4.1). 

2)  The  beam  dimensions  are  much  larger  than  the  width  of 
the  slit.  This  makes  the  one  dimension  analysis  justified. 

3)  The  Debye  length  A^  and  the  ion-electron  mean  free 
32  33 

path  A  .  '  are  both  much  smaller  than  the  size  of  the  slit, 

ei 

Typical  values  are: 

A_  =  4.7x10  6  cm  A  .  =  2.15x10  ^cm  at  N  =  lO^cm  3  and 
D  ei  e 

T  =  4  eV 
e 

A_  =  1.49x10  3 cm  A  .  =  2.15x10  2cra  at  N  =  1017cm  3  and 
D  ei  e 

T  =  40  eV 
e 

4)  The  plasma  can  be  approximated  by  an  inviscid  fluid 

model. 

5)  The  stagnation  of  the  plasma  at  the  plate  surface  is 
consistent  with  our  fluid  theoretical  approach. 

6)  The  steady  state  assumption  is  justified  since  the  beam- 
hole  interaction  time  (<  1  ysec)  is  much  shorter  than  the  lifetime 

of  the  beam  (10  iisec)  and  the  persistence  time  of  the  flow  beyond  the 
hole. 
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7)  Various  photographic  data  analyzing  the  emergence  of  the 
beam  from  the  hole  revealed  that  the  plasma  expands  with  a  measurable 
angle  as  it  comes  out  from  the  opening,  thus  indicating  the  assumed 
compressible  character  of  the  flowing  plasma.  An  incompressible  fluid 
would  completely  fill  in  the  space. 

8)  From  the  temperature  and  velocity  measurements,  one  infers 
that  the  incident  beam  flow  Mach  number  M  =  2  is  raised  to  M  =  4  when 
the  plate  is  present.  These  results  show  that  the  presence  of  the 
plate  changes  the  character  of  the  beam  so  as  to  cool  the  plasma 
considerably  and  at  the  same  time  reduce  its  velocity. 

9)  The  image  converter  camera  pictures  which  provided  an 
integrated  light  exposure  over  the  whole  Mach  wave  pattern,  did  not 
provide  a  check  on  the  theoretical  value  of  the  initial  flow  angle 
0^  as  a  function  of  the  Mach  number  M  {Figure  2.9).  These  data, 
however,  provided  a  means  of  finding  experimentally  the  expanding 
flow  angle  02  provided  the  measurements  are  made  at  the  earliest 
stages  of  the  jet  development.  The  experimental  values  of  ©2  were 
slightly  smaller  than  the  theoretically  predicted  expanding  angles. 
Since  only  visible  light  is  recorded  on  the  Polaroid  film,  the  fast 
cooling  of  the  expanding  plasma  jet  especially  near  the  expansion 
angle  02  will  not  be  recorded  on  the  film,  thus  giving  an  experi¬ 
mental  angle  ©2  less  than  the  one  predicted  by  the  theory. 

Framing  pictures  of  the  beam-hole  interaction  (Fig.  4.19)  taken 
at  latter  stages  in  the  development  of  the  jet  indicated  large  values 
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at  3  ysec 


of  0^  such  as 
02 

02  =  71°  at  5  ysec 

This  is  explained  by  the  fact  that  as  the  time  goes  on,  the 
beam  velocity  after  the  hole  decreases  significantly  as  shown  in 
Fig.  4.16,  while  the  change  in  temperature  is  insignificant  as  illus¬ 
trated  in  Fig.  4.18,  thus  causing  the  Mach  number  to  decrease .  This 
verifies  the  predicted  behavior  of  0^  getting  larger  for  low  Mach 
numbers  as  shown  in  Fig.  2.9. 

10)  The  initial  Mach  number  M  may  be  changed  by  changing  the 
initial  beam  velocity  and  temperature  by  varying  the  gun  voltage  and 
the  back  pressure.  The  relation  between  these  different  parameters 
are  shown  in  Figures  4.3,  4.4,  and  4.5.  In  particular,  the  angles 
6^  and  02  will  thus  be  strongly  affected  by  the  initial  Mach  number 
as  illustrated  in  Fig.  2.9. 

5.2  Conclusion 

As  the  jet  travels  away  from  the  hole,  the  rise  of  the  beam 
temperature  with  distance  {Fig.  4.18)  seems  to  contradict  the  theore¬ 
tical  predictions  of  the  drop  in  temperature  as  shown  in  Figures  2.7 
and  2.8.  This  slight  rise  of  temperature  may  be  attributed  to  weak 
shock  wave  formation  after  the  hole,  which  cannot  be  detected  with 


the  image  converter  camera  due  to  overshadowing  of  the  original  plasma 
light  by  the  ablated  material.  Since  initial  shock  wave  formation 
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was  excluded  in  the  theoretical  considerations,  the  formation  of  shock 
wave  is  explained  as  follows.  At  the  early  stages  of  the  jet  escape, 
the  ratio  between  the  exit  pressure  of  the  jet  and  the  back  pressure 
of  the  outer  medium  is  still  very  high  (“)  causing  the  jet  to  expand 
with  a  given  diverging  angle  (02  <  40°) •  At  later  stages,  due  to  the 
slow  evacuation  of  the  vacuum  chamber  by  the  vacuum  pumps,  the  back 
pressure  in  the  outer  medium  increases  thus  causing  the  exit  to  back 
pressure  ratio  to  become  finite,  leading  to  the  possibility  of 
generating  compression  waves  eventually  leading  to  the  formation  of 
shock  waves  on  the  escaping  jet. 

When  the  effects  of  more  than  one  hole  are  considered,  the 
theoretical  calculations  carried  out  in  Chapter  2  predicted  a  highly 
mixed  plasma  region  in  the  vicinity  of  the  holes  at  low  Mach  numbers 
(Fig.  2.11)  mainly  due  to  the  large  diverging  angle  The  extent 

of  the  incident  plasma  parameters,  i.e.,  density  and  temperature  are 
confined  to  a  triangular  region  (Fig.  2.10),  having  large  cone  angles 
for  low  Mach  numbers.  Therefore,  a  non -uniformity  of  the  plasma 
density  is  expected  in  both  axial  and  azimuthal  directions  in  the 
implosion  chamber.  On  the  contrary,  at  high  Mach  numbers,  very  small 
plasma  mixing  is  expected  to  take  place  within  the  height  of  the 
implosion  chamber.  However,  in  these  cases,  the  incident  plasma 
state  extends  more  into  the  implosion  chamber.  Finally,  a  continuity 
of  the  incident  plasma  beam  state  emerges  out  of  the  holes  forming 
plasma  bars  at  extremely  large  Mach  numbers  due  to  the  presence  of 


the  90  holes. 
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Periodic  azimuthal  non-uniformity  is  therefore,  expected  for 
very  large  Mach  numbers. 

One  recalls  that  in  Chapter  I  the  calculated  minimum  number 
density  required  for  the  necessary  plasma  injection  has  been  deter¬ 
mined  from  the  ideal  case  of  an  infinite  Mach  number.  From  this 
discussion  it  appears  that  an  optimum  Mach  number  can  be  found  that 
is  needed  to  fulfill  the  requirements  of  the  Shiva  experiment,  as 
far  as  the  uniformity  considerations  are  concerned. 

For  the  Mach  number  M  =  4  used  in  the  experiments,  theoretical 

calculations  predict  that  an  initial  plasma  density  N  of  the  order 

e 

of  *v.  2.43xl017cm  3  is  necessary  to  assure  at  least  a  number  density  of 
17  -3 

^  10  cm  or  more  from  each  hole  within  a  slab  of  1  cm  wide  and  2  cm 
long.  Thus,  the  plasma  density  produced  by  the  gun  operating  at  gun 
voltage  't  11  kV,  valve  voltage  'v  8  kV  and  back  pressure  of  'v,  2  psi 
fulfills  the  above  criteria  for  minimum  density  in  a  Shiva  implosion 
chamber . 

Also,  for  Mach  number  4  and  an  opening  of  1  cm,  the  initial 
state  extends  to  a  distance  X^  3.87  cm  falling  therefore  beyond 
the  length  of  the  implosion  chamber  indicating  that  no  reflection  of 
the  expansion  Mach  waves  for  the  axial  wall  symmetry  is  expected. 

One  notices  that,  if  the  separation  between  two  neighboring  holes  is 
1  cm,  the  plasma  mixing  occurs  very  close  to  the  entry  holes 
X2  *\»  0.6  cm. 


With  the  existing  diagnostic  equipment,  due  to  high  level  of 


optical  radiation  by  the  ablated  materials,  the  spatial  variations  of 


the  plasma  density  after  the  hole  were  not  measured. 

Visual  observations  revealed  that  considerable  amounts  of  elec¬ 
trode  material  has  been  ablated  from  the  hole  region  and  transported 
by  the  incident  plasma  beam  into  the  implosion  chamber.  In  fact,  from 
simple  geometrical  considerations,  it  was  found  that  the  ablatant 
material,  injected  into  the  implosion  chamber  in  addition  to  the 
original  He  plasma,  is  evaluated  to  be  2  mgms  (for  Al)  per  shot. 

This  ablation  phenomenon  may  change  the  concept  of  the  plasma  injec¬ 
tion  in  the  implosion  chamber.  With  such  high  ablation  from  high  Z 
materials,  the  ablatants  may  thus  provide  more  than  the  required  mass 
density  for  the  implosion  chamber,  thus  considerably  relaxing  the 
injection  criteria  for  a  purely  gaseous  plasma. 

5. 3  Suggestions  for  Further  Investigations 

In  order  to  fully  understand  the  flow  of  a  plasma  beam  through 
a  hole,  further  study  is  required  both  theoretically  and  experi¬ 
mentally. 

One  of  the  main  investigations  will  be  to  clarify  in  more 
detail  the  formation,  structure  and  dynamics  of  the  beam  ejected 
from  the  coaxial  gun.  For  this,  some  modifications  have  to  be  under¬ 
taken  on  the  gun  in  order  to  facilitate  the  experimental  measurements 
inside  the  barrel  of  the  gun. 

In  the  study  of  the  beam-hole  interaction  itself  a  very  impor¬ 


tant  parameter  which  needs  to  be  measured  experimentally  is  the 
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distribution  of  the  magnetic  field.  The  presence  of  the  magnetic  field 
and  of  electric  currents  flowing  between  the  interaction  plate  and 
the  gun  will  modify  the  fluid  equations  used  in  Chapter  2  as  well  as 
the  dynamics  of  the  beam  formation  in  the  gun. 

In  order  to  visualize  the  Mach  wave  pattern  in  detail,  Schlieren 
34 

photographic  techniques  are  recomnended.  The  eventual  existence 
of  any  shock  wave  at  the  escaping  jet  structure  will  be  apparent  by 
this  technique. 

Also  spatial  and  temporal  distribution  of  the  temperature  and 
density  after  the  opening  are  important  parameters  to  be  monitored 
in  order  to  fully  analyze  the  escaping  jet.  To  avoid  a  single  shot 
type  experiment,  a  sampling  unit  should  be  incorporated  in  the 
diagnostic  equipment  to  evaluate  line  profile  and  continuum  measure¬ 
ments  simultaneously. 

Another  important  further  study  will  be  the  systematic  inves¬ 
tigations  of  the  mechanism  of  the  electrode  material  ablation. 
Specifically,  the  temporal  and  spatial  characteristics  of  the  ablatants 
are  needed.  A  double  holographic  interferometry  technique"^  using 
the  line  wings  of  the  optical  radiation  from  the  ablatant  materials 
is  suggested  for  that  purpose.  For  the  evaluation  of  the  amount  of 
ablatants  deposited  on  the  second  electrode,  a  thin  tilr  interfero¬ 
metry  may  be  used. 

Also,  as  pointed  out  in  Chapter  2,  for  the  plasma  beam-hole 
interaction,  cylindrical  coordinates  should  be  used  in  the  theory 


in  order  to  obtain  more  realistic  solutions  to  the  injection  problem. 
Study  should  be  focussed  on  a  computer  solution  to  carry  out  the 
necessary  iterations  in  the  computations. 
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